1 0 C ACHING OF INTRA-LAYER CALCULATIONS FOR 

RAPID RIGOROUS COUPLED-WAVE ANALYSES 

15 Related Documents 

The present patent application is based on provisional patent application serial 
number 60/178,91 0, filed January 26, 2000, by Xinhui Niu and Nickhil Harshavardhan 
Jakatdar, entitled Cached Coupled Wave Method for Diffraction Grating Profile 
Analysis. 

20 

Background of the Invention 

The present invention relates generally to the measurement of complete profiles 
of one-dimensionally and two-dimensionally repeating, regularly-spaced series of 
structures through the use of an algorithm for caching of intermediate results, and the use 

25 of cached intermediate results to increase the efficiency of calculations. The present 
invention also relates to coupled wave analyses of diffracted electromagnetic signals 
from periodic gratings, particularly where the diffracted electromagnetic signals are 
monitoredby optical metrology tools, such as reflectometer and ellipsometers. More 
particularly, the present invention relates to apparatus and methods for measurement of 

30 complete profiles of one-dimensionally and two-dimensionally repeating, regularly- 
spaced series of structures and to reducing the computation time of coupled wave 
analyses of diffraction from periodic gratings, and still more particularly the present 
invention relates to apparatus and methods for caching and retrieval of intermediate 
computations to reduce the computation time of coupled wave analyses of diffraction 

3 5 from periodic gratings. 
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Diffraction gratings have been used in spectroscopic applications, i.e., 
diffraction applications utilizing multiple wavelengths, such as optical instruments, space 
optics, synchrotron radiation, in the wavelength range from visible to x-rays. 
Furthermore, the past decades have seen the use of diffraction gratings in a wide variety 
5 of nonspectroscopic applications, such as wavelength selectors for tunable lasers, beam- 
sampling elements, and dispersive instruments for multiplexers. 

Advances in computing power are the result of increased speed of semiconductor 
devices. This has been achieved by continually reducing the transistor linewidths, ie. 9 the 
dimensions of transistors. However, as the scale of semiconductor devices decreases, 
10 control of the complete profile of the features, rather than the linewidth alone, is crucial 
for effective chip operation. When the sidewalls are perfectly vertical, knowledge of 
linewidth (and height) is enough to reconstruct the entire line, ie., two parameters define 
the structure. However, due to limits in current fabrication technologies, vertical 
sidewalls are difficult, if not impossible, to achieve, and sloping sidewalls are common. 
15 In addition to the slope of a sidewall, other features that are artifacts of the technology 
which may be present in lines include T-topping (the formation of a "T" shaped profile) 
and footing (the formation of an inverse "T" shaped profile). Capturing such details about 
the profile is important in achieving a better understanding the fabrication technologies. 
In addition to measuring such features, controlling them is also important in this highly 
n 20 competitive marketplace. There are thus increasing efforts to develop and refine run-to- 
run and real-time fabrication control schemes that attempt to make profile measurements 
in-line or in-situ, and use this information to reduce process variability. 

The use of reflectance metrology for the measurement of the thickness of films is 
well-known. In reflectance metrology, an unpolarized or polarized beam of broadband 
25 light is directed towards a sample, and the reflected light is collected. The reflectance can 
either be measured as absolute value, or relative value when normalized to some 
reflectance standards. The reflectance signal is then analyzed to determine the 
thicknesses and optical constants of the film or films. There are numerous examples of 
reflectance metrology. For example, U.S. Pat. No. 5,835,225 given to Thakur et.al. 
30 teaches the use of reflectance metrology to monitor the thickness and refractive indices of 
a film. 
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The use of ellipsometry for the measurement of the thickness of films is also well- 
known (see, for instance, R. M. A. Azzam and N. M. Bashara, "Ellipsometry and 
Polarized Light", North Holland, 1987). When ordinary, i.e., non-polarized, white light 
is sent through a polarizer, it emerges as linearly polarized light with its electric field 
5 vector aligned with an axis of the polarizer. Linearly polarized light can be defined by 
two vectors, i.e., the vectors parallel and perpendicular to the plane of incidence. 
Ellipsometry is based on the change in polarization that occurs when a beam of polarized 
light is reflected from a medium. The change in polarization consists of two parts: a 
phase change and an amplitude change. The change in polarization is different for the 

10 portion of the incident radiation with the electric vector oscillating in the plane of 
incidence, and the portion of the incident radiation with the electric vector oscillating 
perpendicular to the plane of incidence. Ellipsometry measures the results of these two 
changes which are conveniently represented by an angle A, which is the change in phase 
of the reflected beam p from the incident beam; and an angle which is defined as the 

15 arctangent of the amplitude ratio of the incident and reflected beam, Le., 

p = ^ = tanOF)e' (A) , 

where r p is the p-component of the reflectance, and r s is the s-component of the 
reflectance. The angle of incidence and reflection are equal, but opposite in sign, to each 
other and may be chosen for convenience. Since the reflected beam is fixed in position 

20 relative to the incident beam, ellipsometry is an attractive technique for in-situ control of 
processes which take place in a chamber. 

There are numerous examples of the use of spectroscopic ellipsometry. For 
example, U.S. Pat. No. 5,131,752 by Yu et.al. teaches the use of ellipsometry to monitor 
the thickness of a film as it is deposited on a workpiece. The method is however limited 

25 to planar surfaces. U.S. Pat. No. 5,739,909 by Blayo et.al. teaches a method for using 
spectroscopic ellipsometry to measure linewidths by directing an incident beam 131 of 
polarized light at aperiodic structure 100, which has a pitch comprising a feature 121 and 
a space, as is depicted in FIG, 1 . A diffracted beam 132, which leaves the periodic 
structure at an angle 9 which is equal, but opposite in sign, to the angle 6 the incident 

30 beam makes with the periodic structure, is detected and its intensity and polarization are 
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determined at one or more wavelengths. This is then compared with either pre-computed 
libraries of signals or to experimental data, to extract linewidth information. While this is 
a non-destructive test, it does not provide profile information, but yields only a single 
number to characterize the quality of the process. Another method for characterizing 
5 features of a patterned material is disclosed in U.S. Patent No. 5,607,800 by D. H. Ziger. 
According to this method, the intensity, but not the phase, of zeroth-order diffraction is 
monitored for a number of wavelengths, and correlated with features of the patterned 
material. 

While numerous non-destructive techniques have been suggested for linewidth 

10 measurements, such as the scanning electron microscope (SEM) and optical microscope, 
none of them have the ability to provide complete profile information. There exist cross- 
sectional profile metrology tools, such as the atomic force microscope (AFM) and the 
transmission electron microscope, that provide profile information, but at the cost of 
being either prohibitively slow or destructive. Further disadvantages include that these 

1 5 techniques cannot be implemented in-line or in-situ. Finally, there exists scatterometry 
techniques, such as U.S. Pat. No. 5,867,276 by McNeil et.al, which teaches a method for 
measuring profile information. This is accomplished by directing multiple wavelength, 
polarized light onto a periodic structure at a single angle of incidence, and collecting the 
diffracted intensity signal. It is important to note that the incident beam is of a single 

20 planar polarization. This diffracted signal is then compared to a pre-compiled library of 
signals to extract profiles of features. The library is either pre-computed based on 
theoretical calculations or is based purely on experimental signals. This method uses only 
the intensity of the optical signal and has been shown to suffer from non-uniqueness, i.e., 
there exist scenarios where two completely different profiles yield the same intensity 

25 signal, even across a broad range of wavelengths (see, for example, S. Bushman, S. 

Farrer, "Scatterometry Measurements for Process Monitoring of Gate Etch", AEC/APC 
Workshop IX, Sematech, September 20-24 1997). This non-uniqueness reduces 
robustness and accuracy of the results. 

The ability to determine the diffraction characteristics of periodic gratings 
30 with high precision is useful for the refinement of existing applications. Furthermore, the 
accurate determination of the diffraction characteristics of periodic gratings is useful in 
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extending the applications to which diffraction gratings may be applied. However, it is 
well known that modeling of the diffraction of electromagnetic radiation from periodic 
structures is a complex problem that requires sophisticated techniques. Closed analytic 
solutions are restricted to geometries which are so simple that they are of little interest, 
5 and current numerical techniques generally require a prohibitive amount of computation 
time. 

The general problem of the mathematical analysis of electromagnetic diffraction 
from periodic gratings has been addressed using a variety of different types of analysis, 
and several rigorous theories have been developed in the past decades. Methods using 

10 integral formulations of Maxwell's equations were used to obtain numerical results by A. 
R. Neureuther and K. Zaki ("Numerical methods for the analysis of scattering from 
nonplanar periodic structures/' Intn 7 URSI Symposium on Electromagnetic Waves, 
Stresa, Italy, 282-285, 1969) and D. Maystre ("A new general integral theory for 
dielectric coated gratings" 1 Opt Soc. Am., vol. 68, no. 4 5 490-495, April 1978). 

1 5 Methods using differential formulations of Maxwell's equations have also been 
developed by a number of different groups. For instance, an iterative differential 
formulation has been developed by M. Neviere, P. Vincent, R. Petit and M. Cadilhac 
("Systematic study of resonances of holographic thin film couplers, " Optics 
Communications, vol. 9, no. 1, 48-53, Sept. 1973), and the rigorous coupled-wave 

20 analysis method has been developed by M. G. Moharam and T. K. Gaylord ("Rigorous 
Coupled- Wave Analysis of Planar-Grating Diffraction," J. Opt. Soc Am,, vol 71, 81 1- 
818, July 1981). Further work in RCWA formulations has been done by E. B. Grann and 
D. A. Pommet ("Formulation for Stable and Efficient Implementation of the Rigorous 
Coupled-Wave Analysis of Binary Gratings " 1 Opt Soc. Am. A, vol. 12, 1068-1076, 

25 May 1995), and R B. Grann and D. A. Pommet ("Stable Implementation of the Rigorous 
Coupled- Wave Analysis for Surface-Relief Dielectric Gratings: Enhanced Transmittance 
Matrix Approach", J. Opt Soc. Am. A, vol. 12, 1077-1086, May 1995). 

Conceptually, an RCWA computation consists of four steps: 

• The grating is divided into a number of thin, planar layers, and the section of the 
30 ridge within each layer is approximated by a rectangular slab. 
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• Within the grating, Fourier expansions of the electric field, magnetic field, and 
permittivity leads to a system of differential equations for each layer and each 
harmonic order. 

0 Boundary conditions are applied for the electric and magnetic fields at the layer 
5 boundaries to provide a system of equations. 

• Solution of the system of equations provides the diffracted reflectivity from the 
grating for each harmonic order. 

The accuracy of the computation and the time required for the computation depend on the 
number of layers into which the grating is divided and the number of orders used in the 
10 Fourier expansion. 

A number of variations of the mathematical formulation of RCWA have been 
S proposed. For instance, variations of RCWA proposed by P. Lalanne and G. M. Morris 
^ ("Highly Improved Convergence of the Coupled- Wave Method for TM Polarization," J. 
a Opt. Soc. Am. A, 779-784, 1996), L. Li and C. Haggans ("Convergence of the coupled- 
1 5 wave method for metallic lamellar diffraction gratings", J. Opt. Soc. Am. A, 11 84-1 1 89, 
"■3 June, 1993), and G. Granet and B. Guizal ("Efficient Implementation of the Coupled- 
^ Wave Method for Metallic Lamellar Gratings in TM Polarization", J. Opt. Soc. Am. A, 

|t4 1019-1 023 , May, 1 996) differ as whether the Fourier expansions are taken of the 

]fi permittivity or the reciprocal of the permittivity. (According to the lexography of the 
H 20 present specification, all of these variations are considered to be "RCWA.") For a 
specific grating structure, there can be substantial differences in the numerical 
convergence of the different formulations due to differences in the singularity of the 
matrices involved in the calculations, particularly for TM-polarized and conically- 
polarized incident radiation. Therefore, for computational efficiency it is best to select 
25 amongst the different formulations. 

Frequently, the profiles of a large number of periodic gratings must be 
determined. For instance, in determining the ridge profile which produced a measured 
diffraction spectrum in a scatterometry application, thousands or even millions of profiles 
must be generated, the diffraction spectra of the profiles are calculated, and the calculated 
30 diffraction spectra are compared with the measured diffraction spectrum to find the 
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calculated diffraction spectrum which most closely matches the measured diffraction 
spectrum. Further examples of scatterometry applications which require the analysis of 
large numbers of periodic gratings include U.S. Patent Nos. 5,164,790, 5,867,276 and 
5,963,329, and "Specular Spectroscopic Scatterometry in DUV lithography," X. Niu, N. 
Jakatdar, J. Bao and C.J. Spanos, SPIR vol. 3677, pp. 159-168, from thousands to 
millions of diffraction profiles must be analyzed. However, using an accurate method 
such as RCWA, the computation time can be prohibitively long. Thus, there is a need for 
methods and apparatus for rapid and accurate analysis of diffraction data to determine the 
profiles of periodic gratings. 

Additional objects and advantages of the present application will become 
apparent upon review of the Figures, Detailed Description of the Present Invention, and 
appended Claims. 



Summary of the Invention 

The present invention is directed to a method for reducing the computation 
time of rigorous coupled-wave analyses (RCWA) of the diffraction of electromagnetic 
radiation from a periodic grating. RCWA calculations involve the division of the 
periodic grating into layers, with the initial layer being the atmospheric space above the 
grating, the last layer being the substrate below the grating, and the periodic features of 
the grating lying in intermediate layers between the atmospheric space and the substrate. 
A cross-section of the periodic features is discretized into a plurality of stacked 
rectangular sections, and within each layer the permittivity, and the electric and magnetic 
fields of the radiation are formulated as a sum of harmonic components along the 
direction of periodicity of the grating. 

Application of Maxwell's equations provides an intra-layer matrix equation in 
each of the intermediate layers / of the form 



= MK] 
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where Si, y are harmonic amplitudes of an electromagnetic field, z is the perpendicular to 
the periodic grating, and the wave-vector matrix A\ is only dependent on intra-layer 
parameters and incident-radiation parameters. A homogeneous solution of the intra-layer 
matrix equation involves an expansion of the harmonic amplitudes St >y into exponential 
5 functions dependent on eigenvectors and eigenvalues of said wave-vector matrix A\. 

According to the present invention, a layer-property parameter region, an 
incident-radiation parameter region, a layer-property parameter-region sampling, and an 
incident-radiation parameter-region sampling are determined. Also, a maximum 
harmonic order to which the electromagnetic fields are to be computed is determined. 

10 The required permittivity harmonics are calculated for each layer-property value, as 
determined by the layer-property parameter-region sampling of the layer-property 
parameter region. The wave-vector matrix A and its eigenvectors and eigenvalues are 
calculated for each layer-property value and each incident-radiation value, as determined 
by the incident-radiation parameter-region sampling of the incident-radiation parameter 

15 region. The calculated eigenvectors and eigenvalues are stored in a memory for use in 
analysis of the diffraction of incident electromagnetic radiation from the periodic grating. 



Brief Description of the Figures 

Figure 1 shows a section of a diffraction grating labeled with variables used in the 
20 mathematical analysis of the present invention. 

Figure 2 shows a cross-sectional view of a pair of ridges labeled with dimensional 
variables used in the mathematical analysis of the present invention. 

Figure 3 shows a process flow of a TE-polarization rigorous coupled-wave 
analysis. 

25 Figure 4 shows a process flow for a TM-polarization rigorous coupled- wave 

analysis. 

Figure 5 shows a process flow for the pre-computation and caching of calculation 
results dependent on intra-layer and incident-radiation parameters according to the 
method of the present invention. 
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Figure 6 shows a process flow for the use of cached calculation results dependent 
on intra-layer and incident-radiation parameters according to the method of the present 
invention. 

Figure 7 A shows an exemplary ridge profile which is discretized into four stacked 
5 rectangular sections. 

Figure 7B shows an exemplary ridge profile which is discretized into three 
stacked rectangular sections, where the rectangular sections have the same dimensions 
and x-offsets as three of the rectangular section found in the ridge discretization of FIG. 
7A. 

10 Figure 8 shows a computation system for implementation of the computation 

portion of present invention. 

Figure 9 is a flowchart depicted the sequence of events for the implementation of 
the method of the present invention. 

Figure 1 0 shows an apparatus for producing radiation incident on a periodic 
1 5 grating and monitoring the radiation diffracted from the periodic grating. 

Figure 1 1 shows an apparatus for producing radiation incident on a periodic 
grating and monitoring the radiation diffracted from the periodic grating at two angles of 
incidence and diffraction. 

20 

Detailed Description of the Preferred Embodiment 

The method and apparatus of the present invention dramatically improves the 
performance and efficiency for metrology of profiles of one-dimensionally or two- 
dimensionally repeating, regularly-spaced series of structures, by reducing the 
25 computation time required for RCWA computations by pre-processing and caching intra- 
layer information and incident-radiation information, and using the cached intra-layer and 
incident-radiation information for RCWA calculations. 

Section 1 of the present Detailed Description describes apparatus and methods for 
the acquisition of diffraction data from one-dimensionally and two-dimensionally 
30 repeating, regularly-spaced series of structures using optical metrology. 
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Section 2 of the present Detailed Description describes the mathematical 
formalism for RCWA calculations for the diffraction of TE-polarized incident radiation 
from a periodic grating. Definitions of the variables used in the present specification are 
provided, and intra-iayer Fourier-space versions of Maxwell's equations are presented 

5 and solved, producing z-dependent electromagnetic-field harmonic amplitudes, where z is 
the direction normal to the grating. Formulating the electromagnetic-field harmonic 
amplitudes in each layer as exponential expansions produces an eigenequation for a 
wave-vector matrix dependent only on intra-layer parameters and incident-radiation 
parameters. Coefficients and exponents of the exponential harmonic amplitude 

10 expansions are functions of the eigenvalues and eigenvectors of the wave-vector 

matrices. Application of inter-layer boundary conditions produces a boundary-matched 
system matrix equation, and the solution of the boundary-matched system matrix 
equation provides the remaining coefficients of the harmonic amplitude expansions. 

Section 3 of the present Detailed Description describes mathematical formalisms 

1 5 for RCWA calculations of the diffracted reflectivity of TM-polarized incident radiation 
which parallels the exposition of Section 1. 

Section 4 of the present Detailed Description presents a preferred method for the 
solution of the boundary-matched system matrix equation. 

Section 5 of the present Detailed Description describes the method and apparatus 

20 of the present invention. Briefly, the pre-calculation/caching portion of the method of the 
present invention involves: 

• selection of an intra-layer parameter region, an intra-layer parameter sampling, an 
incident-radiation parameter region, and an incident-radiation parameter 
sampling; 

25 • generation of wave- vector matrices for intra-layer parameters spanning the intra- 
layer parameter region, as determined by the intra-layer parameter sampling, and 
incident-radiation parameters spanning the incident-radiation parameter region, as 
determined by the incident-radiation parameter sampling; 

• solution for the eigenvectors and eigenvalues of the wave- vector matrices in the 
3 0 investigative region; and 

• caching of the eigenvectors and eigenvalues of the wave-vector matrices. 
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Briefly, the portion of the method of the present invention for the use of the cached 
computations to calculate the diffracted reflectivity produced by a periodic grating 
includes the steps of; 

• discretization of the profile of a ridge of the periodic grating into layers of 
5 rectangular slabs; 

• retrieval from cache of the eigenvectors and eigenvalues for the wave-vector 
matrix corresponding to each layer of the profile; 

• compilation of the retrieved eigenvectors and eigenvalues for each layer to 
produce a boundary-matched system matrix equation; and 

10 * solution of the boundary-matched system matrix equation to provide the 

diffracted reflectivity. 
Section 6 of the present Detailed Description describes the matching methodology 
^i used in comparing a measured diffraction spectrum to entries in a library of profile- 
. ..1 spectra pairs, in order to determine the profile corresponding to the measured diffraction 
□ 15 spectrum. 

A flowchart for the complete sequence of events for the present invention is 
%l shown in FIG. 9. First, diffraction signals are acquired 905 by an optical setup, using 

d either a reflectometry or an ellipsometry configuration, as described below and depicted 
r*j in FIGS. 10 and 1 1 . Then, numerical aperture of the focusing optics in the optical path of 
VI 20 the reflectometry or ellipsometry configuration are characterized 91 0, in order to 
^ determine the range of incidence angles that need to be taken into account in the 

calculation of the simulated diffraction spectra. Next, a sub-library of profile shapes is 
constructed 915 based on a range and a sampling of intra-layer parameters characterizing 
the profiles. Intra-layer dependent portions of the rigorous coupled-wave calculation are 
25 then pre-computed 920 for the selected range and sampling of intra-layer parameters, and 
cached 925 in a memory. Then, using the pre-computed and cached intra-layer 
dependent portions of the rigorous coupled-wave calculation, diffraction spectra are 
calculated 930 for each profile in the sub-library, to produce a sub-library of spectra. The 
sub-library of profiles is then indexed 935 with the sub-library of spectra to provide a 
30 library of profile-spectra pairs. Finally, a measured diffraction spectrum obtained from a 
physical profile is compared 940 to the spectra in the library, in order to determine the 
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best-matched calculated spectrum. The profile corresponding to the best-matched 
calculated spectrum is the calculated estimate of the physical profile. 

5 1. Acquisition of Diffraction Data 

As shown in FIG. 10, the incident beam (1010) is generated by an excitation head 
(1020) of an ellipsometer (1030). The incident beam (1010) consists of radiation having 
two polarizations to allow measurement of both intensity and phase of the diffracted 
electromagnetic radiation (1011) monitored by an optical detector (1070). An optical 

10 fiber (1025) transports white light from a source (not shown) to the excitation head 

(1020). The white light could either be polarized light or unpolarized light. The incident 
beam (1010) is directed onto a sample (1045) mounted on a'workpiece (1040) such that 
the angle of incidence 8 measured from the normal n of the sample (1045) is between 20 
and 90 degrees, more preferrably between 30 and 80 degrees, still more preferrably 

15 between 40 and 75 degrees, still more preferrably between 50 and 70 degrees. The reason 
the range from 50 and 70 degrees is preferred is because angles in that range have been 
found to generally be most sensitive for the metrology of grating materials typically 
found in semiconductor processing. Most preferrably, the angle of incidence 0 is near the 
Brewster's angle, i.e., the angle such that the sum of the incident and reflected angles is 

20 90 degrees. When the radiation (1010) is incident at the Brewster's angle, the diffracted 
radiation (1011) consists of only a single polarization. The workpiece (1040) can be 
placed on either a chill plate, a hot plate or the developer module (these modules will 
henceforth be referred collectively as the process plate and assigned the reference 
numeral (1080)) on the wafer track, or in the end chamber of an etcher, or in an end- 

25 station or metrology station in a chemical mechanical polishing tool. The reflected beam 
(1011) leaves the workpiece (1040) at an angle 9 which is equal to the angle of incidence 
9, but on the opposite side of the normal n . In the case that the apparatus (1030) of FIG. 
10 is a spectrometer, then the diffracted radiation (1015) is received by the detector 
(1070) which separates the spectra into two polarizations before passing the signal to a 

30 spectrometer (1 090) via an optical fiber (1 027). In the case that the apparatus (1 030) of 
FIG. 10 is a reflectometer, the diffracted radiation (1015) is sent directly to the 
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spectrometer (1090) via the optical fiber (1026). The spectrometer (1090) then sends the 
signal to either a charge coupled device (not shown) or a photomultiplier (not shown) 
which converts the optical signal to an electrical signal which enters a signal processor, 
such as the analysis processor of the present invention depicted in FIG. 8 and described in 

5 detail below. In the case that the apparatus (1030) of FIG. 10 is an ellipsometer, angles 
A and ^ are determined from the numerical signal processing system of FIG, 8. As 
discussed above, the magnitude of the diffracted radiation (1071) provides x ¥ 9 and the 
relative phase of the two polarizations provides A. In the case that the apparatus (1030) 
of FIG. 10 is a reflectometer. the relative intensity is measured. 

10 According to another preferred embodiment of the present invention for a diffraction 

data acquisition system, the broadband ellipsometry or reflectometry apparatus (1030) of 
FIG. 10 is modified to include focusing optics to reduce the spot size to less than the size 
of the test area. Typically, such focusing optics produce illumination regions of 50 \xm x 
50 jim, or smaller, and utilize a pattern recognition module to center the spot in the test 

15 area. 

According to the present invention, a diffraction data acquisition system may utilize 
multiple detectors and multiple excitation heads. For instance, as shown in FIG. 1 1, the 
system (1030') may utilize two excitation heads (1020a) and (1020b), and two detectors 
(1070a) and (1070b). Two broadband incident beams (1 1 10a) and (1 1 10b) are generated 

20 by excitation heads (1 120a) and (1 120b), with each incident beam (1 1 10a) and (1 1 10b) 
including two polarizations of the electromagnetic wave to allow for both intensity and 
phase measurements. As with the single excitation head (1020) apparatus (1030) of FIG. 
10, an optical fiber (1025) transports broadband radiation from a source (not shown) to 
the two excitation heads (1 120a) and (1 120b). A switching mechanism (1099) directs the 

25 broadband radiation alternately to the two optic fiber branches (1025a) and (1 025b) 
which lead to the two excitation heads (1020a) and (1020b), respectively. In the 
preferred embodiment, the switching mechanism (1099) is a wheel with an opening in 
one semicircle, so that rotation of the switching wheel (1099) by 180 degrees produces a 
switching of the optic fiber branch (1025a) or (1025b) to which the radiation is directed. 

30 According to the preferred embodiment of the present invention, the incident beams are 
directed to workpiece (1 140) so that the angles of incidence Gi and 02 are approximately 
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50 and 70 degrees. Other angles of incidence are also contemplated as useful, depending 
on the workpiece, and angles between zero and 90 degrees are possible. Since 
ellipsometric measurements are made for two angles of incidence 9i and 0 2 , the angles 9i 
and 8 2 are chosen to be disparate enough that redundancy of information is not an issue, 

5 while both 8i and 62 are chosen to be close enough to Brewster's angle that the incident 
beam is sensitive to the grating features. The diffracted beams (1015a) and (1015b) leave 
the workpiece at angles Si and 8 2 , which are equal to the angles of incidence 9i and 02, 
but to the opposite side of the normal n . The diffracted beams (101 5a) and (1015b) are 
received by detectors (1070a) and (1070b). In the case that the apparatus (1030 ? ) of FIG, 

10 1 1 is a spectrometer, the diffracted radiation ( 1 0 1 5a) and ( 1 0 1 5b) is received by detectors 
(1070a) and (1070b) which each separates the diffracted radiation (1015a) and (1015b) 
into two polarizations before passing the signal to a spectrometer (1090) via an optical 
fiber (1027). In the case that the apparatus (1030) of FIG. 10 is a reflectometer, then the 
diffracted radiation (1011) is sent directly to the spectrometer (1090) via an optical fiber 

15 (1027). The spectrometer (1090) sends the signal to either a charge coupled device (not 
shown) or a photomultiplier (not shown) which converts the optical signal to an electrical 
signal which enters a signal processor, such as the analysis processor of the present 
invention depicted in FIG, 8 and described in detail below. In the case that the apparatus 
(1030') of FIG, 1 1 is an ellipsometer, angles A and Y as a function of frequency are 

20 determined from the numerical signal processing system of FIG. 8. In the case that the 
apparatus (1030') of FIG. 1 1 is a reflectometer, the relative intensity is measured. It 
should be understood that multiple detectors (1070a) and (1070b) are not necessary to 
perform multiple-angle reflectometry or ellipsometry measurements since the excitation 
head (1020) and detector (1070) may perform a first measurement for an angle of 

25 incidence of Q\ 9 and then be moved to perform a second measurement for an angle of 

incidence of 02. Furthermore, a focusing lenses may be included in the optical path of the 
apparatuses of FIGS. 10 or 11 to produce an area of incidence which small enough to 
meet test-area restrictions. 
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2. Rigorous Coupled- Wave Analysis for TE-Polarized Incident Radiation 

A section of a periodic grating (100) is shown in FIG. 1. The section of the 
grating (100) which is depicted includes three ridges (121) which are shown as having a 
triangular cross-section. It should be noted that the method of the present invention is 
5 applicable to cases where the ridges have shapes which are considerably more complex, 
and even to cases where the categories of "ridges" and "troughs" may be ill-defined. 
According to the lexography of the present specification, the term "ridge" will be used 
for one period of a periodic structure on a substrate. Each ridge (121) of FIG. 1 is 
considered to extend infinitely in the 4-j and -y directions, and an infinite, regularly- 
10 spaced series of such ridges (121) are considered to extend in the -be and -x directions. 
The ridges (121) are atop a deposited film (110), and the film (1 10) is atop a substrate 
(105) which is considered to extend semi-infinitely in the +z direction. The normal 
vector n to the grating is in the -z direction. 

FIG. 1 illustrates the variables associated with a mathematical analysis of a 
15 diffraction grating according to the present invention. In particular: 

• 9 is the angle between the Poynting vector (130) of the incident electromagnetic 
radiation (131) and the normal vector n of the grating (1 00). The Poynting vector 
(130) and the normal vector n define the plane of incidence (140). 

• § is the azimuthal angle of the incident electromagnetic radiation (131), i.e., the 
20 angle between the direction of periodicity of the grating, which in FIG. 1 is along 

the x axis, and the plane of incidence (140). (For ease of presentation, in the 
mathematical analysis of the present specification the azimuthal angle <j> is set to 
zero.) 

• \[/ is the angle between the electric-field vector E of the incident electromagnetic 
25 radiation (131) and the plane of incidence (140), i e., between the electric field 

vector E and its projection E on the plane of incidence (140). When <j>=0 and the 
incident electromagnetic radiation (131) is polarized so that \|/=tt/2, the electric- 
field vector E is perpendicular to the plane of incidence (140) and the magnetic- 
field vector Hlios in the plane of incidence (140), and this is referred to as the TE 



15 



polarization. When <j)=0 and the incident electromagnetic radiation (131) is 
polarized so that \p0, the magnetic-field vector H is perpendicular to the plane of 
incidence (140) and the electric-field vector E lies in the plane of incidence (140), 
and this is referred to as the TM polarization. Any planar polarization is a 
5 combination of in-phase TE and TM polarizations. The method of the present 

invention described below can be applied to any polarization which is a 
superposition of TE and TM polarizations by computing the diffraction of the TE 
and TM components separately and summing them. Furthermore, although the 

'off-axis' (j^O case is more complex because it cannot be separated into TE and 
10 TM components, the present invention is applicable to off-axis incidence 

radiation as well 

* X is the wavelength of the incident electromagnetic radiation (131). 

FIG. 2 shows a cross-sectional view of two ridges (121) of an exemplary 
periodic grating (100) (which will be labeled using the same reference numerals as the 
1 5 grating of FIG. 1 .), illustrating the variables associated with a mathematical description 
of the dimensions of the diffraction grating (100) according to the present invention. In 
particular: 

• L is the number of the layers into which the system is divided. Layers 0 and L are 
considered to be semi-infinite layers. Layer 0 is an "atmospheric" layer (101), 

20 such as vacuum or air, which typically has a refractive index no of unity. Layer L 

is a "substrate" layer (105), which is typically silicon or germanium in 
semiconductor applications. In the case of the exemplary grating 100 of FIG. 2, 
the grating (100) has ten layers with the atmospheric layer (101) being the zeroeth 
layer (125.0), the ridges (121) being in the first through seventh layers (125. 1) 

25 through (125.7), the thin film (110) being the eighth layer (125.8), and the 

substrate (105) being the ninth layer (125.9). (For the mathematical analysis 
described below, the thin-film (1 10) is considered as a periodic portion of the 
ridge (121) with a width d equal to the pitch D) The portion of ridge (121) within 
each intermediate layer (125.1) through (125. (L-l)) is approximated by a thin 
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planar slab (126) having a rectangular cross-section. (Generically or collectively, 
the layers are assigned reference numeral (125), and, depending on context, 
"layers (125)" may be considered to include the atmospheric layer (101) and/or 
the substrate (105).) In general, any geometry of ridges (121) with a cross-section 
which does not consist solely of vertical and horizontal sections can be better 
approximated using a large number of layers (125). 

• Dis the periodicity length or pitch, i e., the length between equivalent points on 
pairs of adjacent ridges (121). 

• di is the width of the rectangular ridge slab (126./) in the /th layer (125./). 

• ti is the thickness of the rectangular ridge slab (126./) in the /th layer (125./) for 
\«{L-\). The thicknesses t\ of the layers (125) are'chosen such that every 
vertical line segment within a layer (125) passes through only a single material. 
For instance, if in FIG. 2 the materials in layers (125.4), (125.5), and (125.6) are 
the same, but different than the materials in layers (125.3) and (125.7), than it 
would be acceptable to combine layers (125.4) and (125.5), or layers (125.5) and 
(125.6), or layers (125.4), (125.5) and (125.6) into a single layer. However, it 
would not be acceptable to combine layers (125.3) and (125.4), or layers (125.6) 
and (125.7) into a single layer. 

• m is the index of refraction of the material in the rectangular ridge slab (1 26) of 
the /th layer (125./). 

In determining the diffraction generated by grating (100), a Fourier space 
version of Maxwell's equations is used. As shown in the calculation process flow 
diagram of FIG. 3, the permittivities si(x) for each layer / are determined or acquired 
(310) (for instance, according to the method described in provisional patent application 
serial number 60/178,540, filed January 26, 2000, entitled Profiler Business Model, by 
the present inventors, and provisional patent application serial number 60/209,424, filed 
June 2, 2000, entitled Profiler Business Model, by the present inventors, both of which 
are incorporated herein by reference), and a one-dimensional Fourier transformation of 
the permittivity e,(x) of each layer / is performed (312) along the direction of periodicity, 
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x , of the periodic grating (100) to provide the harmonic components of the permittivity 
8iu where / is the order of the harmonic component. (In FIGS. 3, 4 5 5 and 6, process steps 
are shown enclosed within ovals or rectangles with rounded corners, and the results of 
calculations are shown enclosed within rectangles with sharp corners. When appropriate 
5 in FIG. 3, equation numbers are used in lieu of, or in addition to, reference numerals.) In 
particular, the real-space permittivity si(x) of the /th layer is related to the permittivity 
harmonics eij of the /th layer by 

s l (x)=^s (tt &xp(j—x). (1.1.1) 
Therefore, via the inverse transform, 
10 s lfi =n 2 r ^+nl(l-^), (1.1.2) 

and for i not equal to zero, 

sm(m — ) 

= ("? -«o) ^~e~^ m , (1.1.3) 

1 5 where n r is the index of refraction of the material in the ridges (12 1 ) in layer /, the index 
of refraction n 0 of the atmospheric layer (101) is typically near unity, and (3 is the x-offset 
of the center of the central rectangular ridge slab (126./) (i.e., the ridge (121) nearest jc=0, 
where generally it is attempted to position the x=0 point at the center of a ridge (121)) 
from the origin. The present specification explicitly addresses periodic gratings where a 

20 single ridge material and a single atmospheric material are found along any line in the re- 
direction. However, as per disclosure document serial number 474051, filed May 15, 
2000, entitled Optical Profilometry for Sub-Micron Periodic Features with Three or More 
Materials in a Layer, by the same inventors, the present invention may be applied to 
gratings having more than one ridge material along a line in the x-direction. 

25 According to the mathematical formulation of the present invention, it is 

convenient to define the (2<H-1) x (2<H-1) Toeplitz-form, permittivity harmonics matrix E\ 
as 
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As will be seen below, to perform a TE-polarization calculation where oth-order 
harmonic components of the electric field E and magnetic field H are used, it is 
5 necessary to use harmonics of the permittivity up to order 2o. 

For the TE polarization, in the atmospheric layer the electric field E is 
formulated (324) as 

£ 0 = expH'Mo ( sin 0x + cos6z) + ^ R, exp(-;(M " K,zi z )) > ( ! - 2 - 1 ) 

i 

where the term on the left of the right-hand side of equation (1 .2.1) is an incoming plane 
1 0 wave at an angle of incidence 9, the term on the right of the right-hand side of equation 
(1.2.1) is a sum of reflected plane waves and R, is the magnitude of the z'th component of 
the reflected wave, and the wave vectors k 0 and (k xi , k 0 , zi ) are given by 

A 

k x =Un,sm{6)-i{^)) , (1-2-3) 

15 and 

_{ k 0 (nj-(kjk 0 ) 2 ) m 24) 

M= l-j* 0 ("MM*o) 2 r ' 

where the value of ^ is chosen from equation (1.2.4), i.e., from the top or the bottom of 
the expression, to provide Re(k 0 , zl ) - Im(% ; ) > 0. This insures that has a positive 
real part, so that energy is conserved. It is easily confirmed that in the atmospheric layer 
20 (101), the reflected wave vector {hi, hzd has a magnitude equal to that of the incoming 
wave vector k o n 0 . The magnetic field #in the atmospheric layer (101) is generated from 
the electric field E by Maxwell's equation (1.3.1) provided below. 
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The x-components k xi of the outgoing wave vectors satisfy the Floquet condition 
(which is also called Bloch's Theorem, see Solid State Physics, N.W. Ashcroft andN.D. 
Mermin, Saunders College, Philadelphia, 1976, pages 133-134) in each of the layers 
(125) containing the periodic ridges (121), and therefore, due to the boundary conditions, 
5 in the atmospheric layer (101) and the substrate layer (105) as well. That is, for a system 
having an n-dimensional periodicity given by 

m = f(r+^m l d l ) , (1-2.5) 
1=1 

where are the basis vectors of the periodic system, and m, takes on positive and 
negative integer values, the Floquet condition requires that the wave vectors k satisfy 



o 



where b, are the reciprocal lattice vectors given by 

(b l -d J ) = S lJ , (1.2.7) 
k 0 is the wave vector of a free-space solution, and 5 tJ is the Kronecker delta function. In 
the case of the layers (125) of the periodic grating (100) of FIGS. 1 and 2 which have the 
1 5 single reciprocal lattice vector b is xlD, thereby providing the relationship of equation 
(1.2.3). 

It may be noted that the formulation given above for the electric field in the 
atmospheric layer (101), although it is an expansion in terms of plane waves, is not 
determined via a Fourier transform of a real-space formulation. Rather, the formulation 

20 is produced (324) a priori based on the Floquet condition and the requirements that both 
the incoming and outgoing radiation have wave vectors of magnitude n 0 k 0 . Similarly, the 
plane wave expansion for the electric field in the substrate layer (105) is produced (324) 
a priori. In the substrate layer (105), the electric field E is formulated (324) as a 
transmitted wave which is a sum of plane waves where the x-components of the wave 

25 vectors {k xi , k 0 , Z i) satisfy the Floquet condition, i e. , 

E^^T^xvHik.x + k^z-^))) , (1.2.8) 
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where 



7*o ("MVW 



,^,^2 fU fu^V2 * (1.2.9) 



where the value of k L)Zi is chosen from equation (1.2.9), ie., from the top or the bottom of 
the expression, to provide Reifa^) - Im(^ zz ) > 0, insuring that energy is conserved. 
5 The plane wave expansions for the electric and magnetic fields in the intermediate 

layers (125.1) through (125.(1-1)) are also produced (334) a priori based on the Floquet 
condition. The electric field E x in the Mi layer is formulated (334) as a plane wave 
expansion along the direction of periodicity, x 9 i.e., 

K = ZM z ) ex P(~^ > 0-2.10) 

1 0 where S t (z) is the z-dependent electric field harmonic amplitude for the /th layer and 
the ith harmonic. Similarly, the magnetic field H t in the /th layer is formulated (334) 
as a plane wave expansion along the direction of periodicity, x , i.e., 

K =-;'(— ) 1/2 I^,.(^)exp(- 7 l ti x) , (1.2.11) 

Mo i 

where U lfXl (z) is the z-dependent magnetic field harmonic amplitude for the Mi layer and 
15 the zth harmonic. 

According to Maxwell's equations, the electric and magnetic fields within a 
layer are related by 

)Vx£, , (1.3.1) 

and 

20 £,=( ~ J )Vxff, . (].3.2) 
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Applying (342) the first Maxwell's equation (1.3.1) to equations (1.2.10) and (1.2.11) 
provides a first relationship between the electric and magnetic field harmonic amplitudes 
Si and Ui of the Zth layer: 
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dS ^ iz) -ku 

ft 0 u I, xi • 

dz 



(1.3.3) 



Similarly, applying (341) the second Maxwell's equation (1.3.2) to equations (1.2.10) and 
(1.2.11), and taking advantage of the relationship 



2nh . 

k xl + = k 

x ' D 



x{,-h) 



(1.3.4) 



which follows from equation (1 .2.3), provides a second relationship between the electric 
and magnetic field harmonic amplitudes Si and U x for the /th layer: 



dz 



K xi 

K k o J 



(1.3.5) 



While equation (1.3.3) only couples harmonic amplitudes of the same order i, equation 
(1 .3.5) couples harmonic amplitudes Si and Ui between harmonic orders. In equation 
(1.3.5), permittivity harmonics s, from order -2o to +2o are required to couple harmonic 
amplitudes Si and U t of orders between -o and +o. 

Combining equations (1.3.3) and (1.3.5) and truncating the calculation to 
order o in the harmonic amplitude S provides (345) a second-order differential matrix 
equation having the form of a wave equation, i.e., 



a 2 s, 



dz'' 



(1.3.6) 



z' = h z, the wave-vector matrix [Aft is defined as 



(1.3.7) 
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where [K x ] is a diagonal matrix with the (i,f) element being equal to (k xi /k 0 ), the 
permittivity harmonics matrix [Eft is defined above in equation (1 .1 .4), and [ ] and 
[cfSiy/dz* 2 ] are column vectors with indices i running from -o to +o, i.e., 



^l,y,0 
$l>y,o 



(1.3.8) 



5 By writing (350) the homogeneous solution of equation (1 3.6) as an 

expansion in pairs of exponentials, i.e., 



2o+\ \ 

S i,yX z ) = X l *W d '.« ex P(-*o^) + c2 />jn expC^^Cz-f;))] , (1.3.9) 



its functional form is maintained upon second-order differentiation by z \ thereby taking 
the form of an eigenequation. Solution (347) of the eigenequation 



(1.3.10) 



provides (348) a diagonal eigenvalue matrix fa] formed from the eigenvalues y tTn of the 
wave-vector matrix [Ai\ 9 and an eigenvector matrix [W{] of entries m//,^, where wi^m is the 
rth entry of the mth eigenvector of [Aft. A diagonal root-eigenvalue matrix [Qj\ is defined 
to be diagonal entries qi tl which are the positive real portion of the square roots of the 
eigenvalues x^. The constants cl and c2 are, as yet, undetermined. 

By applying equation (1.3.3) to equation (1.3.9) it is found that 



2(3+1 - 

U h A z ) =2 v /^r c1 ^ ex P("*o?/,»^) + c 2 /fBI exp(* 0 $ /()n (z-f,))] (1.3.11) 



where v/ t ^ m = qi im Wti t m- The matrix to be used below, is composed of entries v^ m . 

The constants cl and c2 in the homogeneous solutions of equations (1.3.9) 
20 and (1 .3 . 1 1) are determined by applying (355) the requirement that the tangential electric 
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and magnetic fields be continuous at the boundary between each pair of adjacent layers 
(125.0/(125.(7+1)). At the boundary between the atmospheric layer (101) and the first 
layer (125.1), continuity of the electric field E y and the magnetic field H x requires 



I 



R = 



W x W x X x Jc\ 



(1.4.1) 



where Yq is a diagonal matrix with entries (ko iZi / k 0 ), Xi is a diagonal layer-translation 
matrix with elements exp(-fe q^ m tj), R is a vector consisting of entries from R. 0 to R+ 0 
and cl] and c2\ are vectors consisting of entries from cl i } o and clj^h and c2j j0 and 
c2j ! 2o+h respectively. The top half of matrix equation (1 .4. 1) provides matching of the 
electric field E y across the boundary of the atmospheric layer (125.0) and the first layer 
(125.1), the bottom half of matrix equation (1.4.1) provides matching of the magnetic 
field H x across the layer boundary (125.0)/(125.1), the vector on the far left is the 
contribution from the incoming radiation (131) in the atmospheric layer (101), the second 
vector on the left is the contribution from the reflected radiation (132) in the atmospheric 
layer (101), and the portion on the right represents the fields E y and H x in the first layer 
(125.1). 

At the boundary between adjacent intermediate layers (125./) and (125.(/+1)), 
continuity of the electric field E y and the magnetic field H x requires 



cl 



l-l 



cl t _ 



i_ 



W l W l X l ' 

Yi -v,Xi 



cl/ 
c2j 



(1.4.2) 



where the top and bottom halves of the vector equation provide matching of the electric 
20 field Ey and the magnetic field H x , respectively, across the l-l/l layer boundary. 

At the boundary between the (I-l)th layer 125.(L-1) and the substrate layer 
(105), continuity of the electric field E y and the magnetic field H x requires 











I 










J Y L. 



(1.4.3) 
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where, as above, the top and bottom halves of the vector equation provides matching of 
the electric field E y and the magnetic field H x , respectively. In contrast with equation 
(1 .4.1), there is only a single term on the right since there is no incident radiation in the 
substrate (105). 

5 Matrix equation (1.4.1), matrix equation (1 A3), and the (LA) matrix 

equations (1 .4.2) can be combined (360) to provide a boundary-matched system matrix 
equation 



" -/ 
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w x x x 
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jS tQ n Q cos(0) 


0 - 


w x x x 


-w { 


w 2 


W 2 X 2 


0 


0 •■• 


c2 x 




0 


0 - 


v x x x 






-V 7 X 2 


0 


0 ••• 






0 


0 


















0 


0 






-w^x^ 
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T 




0 



(1.4.4) 

10 and this boundary-matched system matrix equation (1 .4.4) may be solved (365) to 

provide the reflectivity R { for each harmonic order i, (Alternatively, the partial-solution 
approach described in "Stable Implementation of the Rigorous Coupled- Wave Analysis 
for Surface-Relief Dielectric Gratings: Enhanced Trausmittance Matrix Approach", E. B. 
Grann and D. A. Pommet, J. Opt Soc. Am. A, vol. 12, 1077-1086, May 1995, can be 

15 applied to calculate either the diffracted reflectivity R or the diffracted transmittance T) 

3. Rigorous Coupled- Wave Analysis for the TM Polarization 

The method (400) of calculation for the diffracted reflectivity of TM-polarized 
incident electromagnetic radiation (131) shown in FIG. 4 parallels that (300) described 
20 above and shown in FIG. 3 for the diffracted reflectivity of TE-polarized incident 

electromagnetic radiation (131). The variables describing the geometry of the grating 
(100) and the geometry of the incident radiation (13 1) are as shown in FIGS. 1 and 2. 
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However, for TM-polarized incident radiation (13 1) the electric field vector £ is in the 
plane of incidence (140), and the magnetic field vector His perpendicular to the plane of 
incidence (140). (The similarity in the TE- and TM-polarization RCWA calculations and 
the application of the present invention motivates the use of the term 'electromagnetic 
5 field' in the present specification to refer generically to either or both the electric field 
and/or the magnetic field of the electromagnetic radiation.) 

As above, once the permittivity sfe) is determined or acquired (410), the 
permittivity harmonics Su are determined (412) using Fourier transforms according to 
equations (1.1.2) and (1.1.3), and the permittivity harmonics matrix Ei is assembled as 
10 per equation (1.1.4). In the case of TM-polarized incident radiation (131), it has been 
found that the accuracy of the calculation may be improved by formulating the 
calculations using inverse-permittivity harmonics since this will involve the inversion 
of matrices which are less singular. In particular, the one-dimensional Fourier expansion 
(412) for the inverse of the permittivity s\(x) of the /th layer is given by 

15 Z exp(y^x) . (2.1.1) 

Therefore, via the inverse Fourier transform this provides 



1 d, 1 d, 

*/.o=- T^ + — (2-1-2) 

nt D n n D 



and for h not equal to zero, 



- - sin(;r/z — ) 

20 x i h = (-y — j) -£- e -"*» D , (2.1.3) 

n r »o 7th 

where p is the x-offset of the center of the rectangular ridge slab (126./) from the origin. 
The inverse-permittivity harmonics matrix Pi is defined as 
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where 2o is the maximum harmonic order of the inverse permittivity n I>h used in the 
calculation. As with the case of the TE polarization (300), for electromagnetic fields E 
and if calculated to order o it is necessary to use harmonic components of the 
permittivity si t h and inverse permittivity %\ } h to order 2o. 

In the atmospheric layer the magnetic field H i? formulated (424) a priori as 
a plane wave incoming at an angle of incidence 9, and a reflected wave which is a sum of 
plane waves having wave vectors (k xl) ko iZl ) satisfmg the Floquet condition, equation 
(1.2.6). In particular, 

H 0j =exp(-/Vo( sinflc + cos& ) + Z^ ex P(-X***-*o,**)) > ( 22A ) 

where the term on the left of the right-hand side of the equation is the incoming plane 
wave, and R t is the magnitude of the zth component of the reflected wave. The wave 
vectors k 0 and (k xl , ko,zt) are S iven b Y equations (1.2.2), (1.2.3), and (1.2.4) above, and the 
magnetic field H in the atmospheric layer (101) is generated from the electric field £by 
Maxwell's equation (1.3.2). In the substrate layer (105) the magnetic field H is 
formulated (424) as a transmitted wave which is a sum of plane waves where the wave 
vectors (k xh k 0fZ i) satisfy the Floquet condition, equation (1.2.6), ie. 9 

i M 

where k LiZi is defined in equation (1.2.9), Again based on the Floquet condition, the 
magnetic field H { in the fth layer is formulated (434) as a plane wave expansion along 

the direction of periodicity, x , i.e., 
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where U, (z) is the z-dependent magnetic field harmonic amplitude for the fth layer and 
the ith harmonic. Similarly, the electric field in the fth layer is formulated (434) as a 
plane wave expansion along the direction of periodicity, i.e., 

E l:X = j(-^-) m (*) exp(-j^x) , (2.2.4) 

where xi {z) is the z-dependent electric field harmonic amplitude for the Mi layer and 
the fth harmonic. 

Substituting equations (2.2.3) and (2.2.4) into Maxwell's equation (1.3.2) 
provides (441) a first relationship between the electric and magnetic field harmonic 
amplitudes Si and Ui for the fth layer: 

^ElA= [El ][S, x ,]. (2.3.1) 
oz 

Similarly, substituting (2.2.3) and (2.2.4) into Maxwell's equation (1 .3.1) provides (442) 
a second relationship between the electric and magnetic field harmonic amplitudes Si and 
Ui for the fth layer: 



flf =([K x ][P l ][K x ]-U])[U l J. (2.3.2) 
oz 



where, as above, K x is a diagonal matrix with the (i, z) element being equal to (k xi /k 0 ). In 
contrast with equations (1.3.3) and (1.3.5) from the TE-polarization calculation, non- 
diagonal matrices in both equation (2.3.1) and equation (2.3.2) couple harmonic 
amplitudes Si and Ui between harmonic orders. 

Combining equations (2.3.1) and (2.3.2) provides a second-order differential 
wave equation 
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8z< 2 



= {[E l ](lK x ][P ! ][K x ]-Ui)}[uJ , (2-3.3) 



where [Ui, y ] and [&U\Jdz a ] are column vectors with indices running from -o to +o, and 
the permittivity harmonics [Eft is defined above in equation (1.1.7), and z'=koz. The 
wave-vector matrix [Aft for equation (2.3.3) is defined as 

5 Wl=[^]taWHfl)- (2-3-4) 

If an infinite number of harmonics could be used, then the inverse of the permittivity 
harmonics matrix [Eft would be equal to the inverse-permittivity harmonics matrix [Pft, 
and vice versa, i.e., [Eft 1 = [Pi], and [Pft' 1 = [Eft. However, the equality does not hold 
when a finite number o of harmonics is used, and for finite 'o the singularity of the 
10 matrices [Eft' 1 and [Pft, and the singularity of the matrices [P,]" 1 and [Eft, will generally 
differ. In fact, it has been found that the accuracy of RCWA calculations will vary 
depending on whether the wave-vector matrix [Aft is defined as in equation (2.3.4), or 

[A^Wr^^tfA-V}), (2-3.5) 



or 



15 [A^&jilK^nKJ-V]). (2.3.6) 

It should also be understood that although the case where 

[A]=^r 1 ([^]^][^J-m) (2-3.6') 

does not typically provide convergence which is as good as the formulations of equation 

(2.3.5) and (2.3.6), the present invention may also be applied to the formulation of 
20 .equation (2.3.6'). 

Regardless of which of the three formulations, equations (2.3.4), (2.3.5) or 

(2.3.6) , for the wave-vector matrix [Aft is used, the solution of equation (2.3.3) is 
performed by writing (450) the homogeneous solution for the magnetic field harmonic 
amplitude £/?as an expansion in pairs of exponentials, i.e., 
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JJ iA z ) = IX,»<K m ex P(-^o^,m z ) + c2 />m expC^^Cz -/,))] . (2.3.7) 

m=l 

since its functional form is maintained upon second-order differentiation by z\ and 
equation (2.3.3) becomes an eigenequation. Solution (447) of the eigenequation 

provides (448) an eigenvector matrix [Wi\ formed from the eigenvectors of the wave- 
vector matrix [A{\, and a diagonal eigenvalue matrix [xj formed from the eigenvalues t/,/ 
of the wave-vector matrix [Aj\. A diagonal root-eigenvalue matrix [Qi\ is formed of 
diagonal entries gu which are the positive real portion of the square roots of the 
eigenvalues The constants c\ and c2 of equation (2.3.7) are, as yet, undetermined. 

By applying equation (1.3.3) to equation (2.3.5) it is found that 

2o+l j. -J 

Si* 00 = Z v '.'.» L~ c1 /.™ ex P(-^o?;, m z ) + c2 i, m exp(* 0 ?^ (* - f, ))J (2.3.9) 

m=l 

where the vectors vij form a matrix [V{\ defined as 

[V] = [£]"* [JF] [g] when [4] is defined as in equation (2.3.4) , (2.3.10) 

[V\ = [P][W\[Q\ when " " ( 2 - 3 - 5 )' ( 23 - n ) 

[V\ = [E\' l W\[Q\ when " " (2.3.6). (2.3.12) 

The formulation of equations (2.3.5) and (2.3.11) typically has improved convergence 
performance (see P. Lalanne and G.M. Morris, "Highly Improved Convergence of the 
Coupled-Wave Method for TM Polarization", J. Opt. Soc. Am. A, 779-784, 1996; and L. 
Li and C. Haggans, "Convergence of the coupled-wave method for metallic lamellar 
diffraction gratings", J. Opt. Soc. Am. A, 1 184-1 189, June 1993) relative to the 
formulation of equations (2.3.4) and (2.3.1 1) (see M.G. Moharam and T.K. Gaylord, 
"Rigorous Coupled-Wave Analysis of Planar-Grating Diffraction", J. Opt. Soc. Am., vol. 
71, 811-818, July 1981). 
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The constants c\ and c2 in the homogeneous solutions of equations (2.3.7) 
and (2.3.9) are determined by applying (455) the requirement that the tangential electric 
and tangential magnetic fields be continuous at the boundary between each pair of 
adjacent layers (125i)/(125.(/+l)), when the materials in each layer non-conductive. 

5 (The calculation of the present specification is straightforwardly modified to 

circumstances involving conductive materials, and the application of the method of the 
present invention to periodic gratings which include conductive materials is considered to 
be within the scope of the present invention. At the boundary between the atmospheric 
layer (101) and the first layer (125.1), continuity of the magnetic field H y and the electric 

10 field E x requires 
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where Z 0 is a diagonal matrix with entries (k 0iZI /n 0 2 k 0 ), X x is a diagonal matrix with 
elements exp(-& 0 qi, m ti), the top half of the vector equation provides matching of the 
magnetic field H y across the layer boundary, the bottom half of the vector equation 
1 5 provides matching of the electric field E x across the layer boundary, the vector on the far 
left is the contribution from the incoming radiation (131) in the atmospheric layer (101), 
the second vector on the left is the contribution from the reflected radiation (132) in the 
atmospheric layer (101), and the portion on the right represents the fields H y and E x in the 
first layer (125.1). 

20 At the boundary between adjacent intermediate layers (125 1) and (125.(/+1)), 

continuity of the electric field E y and the magnetic field H x requires 
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where the top and bottom halves of the vector equation provides matching of the 
magnetic field H y and the electric field E x , respectively, across the layer boundary. 

25 At the boundary between the (Z-l)th layer (125. (L-l)) and the substrate layer 

(105), continuity of the electric field E y and the magnetic field H x requires 
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where, as above, the top and bottom halves of the vector equation provides matching of 
the magnetic field H y and the electric field E x , respectively. In contrast with equation 
(2.4.1), there is only a single term on the right in equation (2.4.3) since there is no 
incident radiation in the substrate (105). 

Matrix equation (2.4.1), matrix equation (2.4.3), and the (14) matrix 
equations (2.4.2) can be combined (460) to provide a boundary-matched system matrix 
equation 
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10 (2.4.4) 

and the boundary-matched system matrix equation (2.4.4) may be solved (465) to provide 
the reflectivity R for each harmonic order I (Alternatively, the partial-solution approach 
described in "Stable Implementation of the Rigorous Coupled- Wave Analysis for 
Surface-Relief Dielectric Gratings; Enhanced Transmittance Matrix Approach", E, B. 
15 Grann and D. A. Pommet, 1 Opt Soc. Am. A, vol. 12, 10774086, May 1995, can be 
applied to calculate either the diffracted reflectivity R or the diffracted transmittance T) 



4. Solving for the Diffracted Reflectivity 

The matrix on the left in boundary-matched system matrix equations (L4.4) 
20 and (2.4.4) is a square non-Hermetian complex matrix which is sparse (Le. 9 most of its 
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entries are zero), and is of constant block construction (i.e., it is an array of sub-matrices 
of uniform size). According to the preferred embodiment of the present invention, and as 
is well-known in the art of the solution of matrix equations, the matrix is stored using the 
constant block compressed sparse row data structure (BSR) method (see S. Carney, M. 
5 Heroux, G. Li, R. Pozo, K. Remington and K. Wu, "A Revised Proposal for a Sparse 
BLAS Toolkit," http://www.netlib.org , 1996). In particular, for a matrix composed of a 
square array of square sub-matrices, the BSR method uses five descriptors: 

BJLDA is the dimension of the array of sub-matrices; 

0 is the dimension of the sub-matrices; 

VAL is a vector of the non-zero sub-matrices starting from the leftmost non-zero 
matrix in the top row (assuming that there is a non-zero matrix in the top row), 
and continuing on from left to right, and top to bottom, to the rightmost non-zero 
matrix in the bottom row (assuming that there is a non-zero matrix in the bottom 
row). 

COLJND is a vector of the column indices of the sub-matrices in the VAL 
vector; and 

ROW__PTR is a vector of pointers to those sub-matrices in VAL which are the 
first non-zero sub-matrices in each row. 

For example, for the left-hand matrix of equation (1 .4.4), B__LDA has a value of 2L, 0 
20 has a value of 2o+l, the entries of VAL are (-7, W h WjX h jY 0t V h -VjX h -WjX h -W }) W 2} 
W 2 X 2) -ViX h V h V 2> ....), the entries of COLJND are (1, 2, 3, 1, 2 5 3, 2, 3, 4, 5, 2, 3, 4, 
5, . . .), and the entries of RO W__PTR are (1 , 4, 7, 1 1 , . . .). 

According to the preferred embodiment of the present invention, and as is 
well-known in the art of the solution of matrix equations, the squareness and sparseness 
25 of the left-hand matrices of equations (1 .4.4) and (2.4.4) are used to advantage by solving 
equations (1.4.4) and (2.4.4) using the Blocked Gaussian Elimination (BGE) algorithm. 
The BGE algorithm is derived from the standard Gaussian Elimination algorithm (see, for 
example, Numerical Recipes, W.H. Press, B.P. Flannery, S.A. Teukolsky, and W.T. 
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Vetterling, Cambridge University Press, Cambridge, 1986, pp. 29-38) by the substitution 
of sub-matrices for scalars. According to the Gaussian Elimination method, the left-hand 
matrix of equation (1.4.4) or (2.4.4) is decomposed into the product of alower triangular 
matrix [L], and an upper triangular matrix [U], to provide an equation of the form 
5 [L][U][x] = [b], (3.1.1) 

and then the two triangular systems [U][x] = [y]anA[L][y] = [b]are solved to 
obtain the solution [ x ] = [ U] A [LJ l [b], where, as per equations (1.4.4) and (2.4.4), 
[ x ] includes the diffracted reflectivity R. 

10 

5. Caching of Permittivity Harmonics and Eigensolutions 

As presented above, the calculation of the diffraction of incident TE-polarized 
or TM-polarized incident radiation (131) from a periodic grating involves the generation 
of a boundary-matched system matrix equation (1 .4.4) or (2.4.4), respectively, and its 

1 5 solution. In understanding the advantages of the present invention it is important to 
appreciate that the most computationally expensive portion of the processes of FIGS. 3 
and 4 is the solution (347) and (447) for the eigenvectors w hi and eigenvalues x hl of wave- 
vector matrix [A,] of equation (1.3.7), (2.3.4), (2.3.5) or (2.3.6). The accuracy of the 
calculation of the eigenvectors and eigenvalues x u is dependent on the number of 

20 orders o utilized. As the number of orders o is increased, the computation time for 
solving the eigensystem increases exponentially with o. When performed in a typical 
computing environment with o=9 harmonic orders, the calculation of the eigenvectors 
and eigenvalues can take more than 85% of the total computation time. 

The method of the present invention is implemented on a computer system 
25 (800) which in its simplest form consists of information input/output (I/O) equipment 
(805), which is interfaced to a computer (810) which includes a central processing unit 
(CPU) (815) and a memory (820). The I/O equipment (805) will typically include a 
keyboard (802) and mouse (804) for the input of information, and a display (801) and 
printer (803) for the output of information. Many variations of this simple computer 
30 system (800) are to be considered within the scope of the present invention, including 
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systems with multiple I/O devices, multiple processors within a single computer, multiple 
computers connected by Internet linkages, multiple computers connected by local area 
networks, devices that incorporate CPU (815), memory (820), and I/O equipment (805) 
into a single unit, etc. For instance, the method of the present invention may be applied 
5 to any of the systems described in provisional patent application serial number 

60/178,540, filed January 26, 2000, entitled Profiler Business Model, by the present 
inventors, and provisional patent application serial number 60/209,424, filed June 2, 
2000, entitled Profiler Business Model, by the present inventors, both of which are 
incorporated herein by reference. 

1 o According to the method and apparatus of the present invention, portions of 

the analysis of FIG. 3 are pre-computed and cached, thereby reducing the computation 
time required to calculate the diffracted reflectivity produced by a periodic grating. 
Briefly, the pre-computation and caching portion of the present invention consists of: 

• pre-computation and caching (i. e. , storage in a look-up table) of the permittivity 
1 5 s^x), the harmonic components s^i of the permittivity e^x) and the permittivity 

harmonics matrix and/or the inverse-permittivity harmonics and the 
inverse-permittivity harmonics matrix [P„] for a sampling region {u} of layer- 
property values; 

• pre-computation and caching of the wave-vector matrix [A^] for the sampling 
20 region { \x} of layer-property values and a sampling region {k} of incident- 
radiation values; and 

• pre-computation and caching of eigenvectors and eigenvalues of the 
wave-vector matrix [A^ K ] to form an eigenvector matrix [W^], a root-eigenvalue 
matrix [Q&k], and a compound matrix [V^ K ], respectively, for a master sampling 

25 region {jj,, k} formed from the combination of the layer-property sampling region 

{u} and the incident-radiation sampling region {k}; 
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Briefly, the use of the master sampling region {\x 9 ic}of pre-computed and cached 
eigenvector matrices [W^ K ], root-eigenvalue matrices [Q^ K ], and product matrices [V^ K ] 
to calculate the diffraction spectrum from a periodic grating consists of the steps of: 

• construction of matrix equation (1 .4.4) or (2.4.4) by retrieval of cached 

5 eigenvector matrices [W^ K ], root-eigenvalue matrices [Q^ K ], and product matrices 

[F^ K ] from the master sampling region {p,, k} corresponding to the layers 125 of 
the grating (100) under consideration; and 

• solution of the matrix equation (1 .4.4) or (2.4.4) to determine the diffracted 
reflectivity R : for each harmonic order i. 

1 0 The method of the present invention is illustrated by consideration of the 

exemplary ridge profiles (701) and (751) shown in cross-section in FIGS. 7 A and 7B, 
respectively. The profile (701) of FIG. 7 A is approximated by four slabs (711), (712), 
(713) and (714) of rectangular cross-section. Similarly, the profile (751) of FIG. 7B is 
approximated by three slabs (761), (762), and (763) of rectangular cross-section. The 

15 two exemplary ridge profiles (701) and (751) are each part of an exemplary periodic 
grating (other ridges not shown) which have the same grating period Z), angle 9 of 
incidence of the radiation (131), and radiation wavelength \. Furthermore, slabs (713) 
and (761) have the same ridge slab width d, x-offset p, and Index of refraction n r> and the 
index of refraction no of the atmospheric material between the ridges (701) and (751) is 

20 the same. Similarly, slabs (711) and (762) have the same ridge slab width d, x-offset p, 
and index of refraction n n and slabs (714) and (763) have the same ridge slab width d, x- 
offset p ? and index of refraction n r . However, it should be noted that slabs (714) and 
(763) do not have the same thicknesses t 9 nor do slabs (713) and (761) or slabs (711) and 
(762) have the same thicknesses t. It is important to note that thickness t is not a 

25 parameter upon which the wave-vector matrix [A] is dependent, although thickness t does 
describe an intra-layer property. It should also be noted that the present invention may be 
implemented with ridges mounted directly on a substrate, or ridges mounted on films 
deposited on a substrate, since a film can be considered to be a ridge having a width d 
equal to the pitch D. 
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In performing an RCWA calculation for the diffracted reflectivity from 
grating composed of profiles (701), the eigenvector matrices [W] 9 the root-eigenvalue 
matrices [Q] 9 and the compound eigensystem matrices [F] are computed for rectangular 
slabs (711), (712), (713), and (714). According to the present invention it is noted that 
5 the eigenvector matrices [W], the root-eigenvalue matrices [Q], and the compound 

eigensystem matrices [F] for slabs (761), (762) and (763) are the same as the eigenvector 
matrices [W] 9 the root-eigenvalue matrices [Q] 9 and the compound eigensystem matrices 
[F] for slabs (713), (711) and (714), respectively, since the wave- vector matrices [A] are 
the same for slabs (71 1) and (762), (713) and (761), and (714) and (763). Therefore, 

10 caching and retrieval of the eigensystem matrices [W] 9 [Q\ and [V] for slabs (713), (71 1) 
and (714) would prevent the need for recalculation of eigensystem matrices [W\, [Q] 9 and 
[V] for slabs (761), (762) and (763), and reduce the computation time. More broadly, the 
pre-calculation and caching of eigensystem matrices [W], [Q] 9 and [F] for useful ranges 
and samplings of intra-layer parameters and incident-radiation parameters can reduce the 

1 5 computation time necessary to perform RCWA calculations. 

As can be seen from equations (1.1.2), (1.1.3), (2.1.2) and (2.13), the 
permittivity harmonics and the inverse permittivity harmonics n ltt are only dependent 
on the intra-layer parameters: the index of refraction of the ridges w^the index of 
refraction of the atmospheric material no 9 the pitch D, the ridge slab width d 9 and the x- 

20 offset p. As shown in the flowchart of FIG. 5, in one exemplary embodiment, the system 
(600) of the present invention begins with the determination (605) of the ranges n ri mn to 
ftr,max 9 no,mm to 7i0 jBMX , D mm to D maX9 d mm to d mQX9 and $ mm to $ maX9 and increments 8n r? Srio, 
Si), 8<tf, and 5p for the layer-property parameters, Le. 9 the index of refraction of the ridges 
n rt the index of refraction of the atmospheric material n 09 the pitch D 9 the ridge slab width 

25 d 9 the x-offset P, as well as the determination (605) of the maximum harmonic order o, 
This information is forwarded from an I/O device (805) to the CPU (815). Typically, 
when applied to periodic gratings produced by semiconductor fabrication techniques, the 
ranges n r>min to n r>ma x> no.mm to n o,max, D mm to D max , d mm to d maX9 and p m/n to $ max are 
determined based on knowledge and expectations regarding the fabrication materials, the 

30 fabrication process parameters, and other measurements taken of the periodic grating 
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(100) or related structures. Similarly, when matching calculated diffraction spectra to a 
measured diffraction spectrum to determine the dimensions of the periodic grating that 
created the measured diffraction spectrum, the increments hn r> §no, SD 5 5d, and 8(3, and 
maximum harmonic order o 3 are chosen based on the resolution to which the layer- 
5 property parameters w r , no, D, d and (3 are to be determined. The layer-property parameter 
ranges n r)min to n r , max , no )min to no )max , D mtn to D max , 4« to d max , and § mn to $ max , and 
increments bn r> S#0, SD, and bd, and 8p define a five-dimensional layer-property caching 
grid {p,}. More specifically, the caching grid {p} consists of layer-property points with 
the TV coordinates being {n rMm n m \n+§n Ti n rimin +2§n rt n rtmesr 28n r9 n r)max -hn r) n rtmax } 9 the 

10 no coordinates being {no, } ? the 
D coordinates being {D mirb D mtn J rhD ) D mm +2hD 9 . . . ? D mwr 2hD t D max -hD, D max } , the d 
coordinates being {d mfu d min J rSd t d min J r2Sd i d max -2Sd i d ma x-Sd r d max }, and the p 
coordinates being {p miM , p mz? ,+<5p, $mm+2d$, Pm^r^, p ma x}« In other words, 

the layer-property caching grid {ji} is defined as a union of five-dimensional coordinates 

15 as follows: 

{/'} = U (n rMn + f , ^ min + jSn, , ^ + *<5D, J mm + /<?J ? £ mn + mSjB) , (4. L 1 ) 

where i, j, k, I and m are integers with value ranges of 

O^K^-^J/^,, (4.1.2a) 

0<j<(n^~n aMa )/dn 0 , (4.1.2b) 

20 O^CD^-I^)/,©, (4.1.2c) 

Q<l<{d^-d^ n )!5d, (4.1.2d) 

and 

OZmZi^-fi^ySfi. (4.1.2.e) 
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It should be noted that the variable / in equations (4. 1 . 1) and (4. 1 .2d) is not to be 
confused with the layer number / used in many of the equations above. Furthermore, it 
may be noted that the layer subscript, /, is not used in describing the layer-property 
parameters n r , no, A 4 P used in the layer-property caching grid {\x} because each 

5 particular point jjj in the layer-property caching grid {ji} may correspond to none, one, 
more than one, or even all of the layers of a particular periodic grating (100). It should 
also be understood that the layer-property parameter region need not be a hyper- 
rectangle, and the layer-property parameter region need not be sampled using a grid. For 
instance, the sampling of the layer-property parameter region may be performed using a 

10 stochastic sampling method. Furthermore, the sampling density of the layer-property 
parameter region need not be uniform. For instance, the sampling density (i.e., the 
sampling resolution) may decrease near the boundaries of the layer-property parameter 
region if layers (125) described by layer properties near the boundaries are less likely to 
occur. 

15 As shown in FIG. 5, for each point juj in the layer-property caching grid {\x} 

the "required" permittivity harmonics s l are calculated 410 by CPU (815) and cached 

(415) in memory (820), and the "required" permittivity harmonics matrices [£]are 

compiled from the cached required permittivity harmonics e l and cached (415') in 

memory (820). For RCWA analyses of TE-polarized incident radiation (131), or RCWA 
20 analyses of TM-polarized incident radiation (131) according to the formulation of 

equations (2.3.6) and (2.3.12), the required permittivity harmonics s i are the permittivity 
harmonics s t calculated (410) according to equations (1.1.2) and (1.1.3), and the required 
permittivity harmonics matrix \e] is the permittivity harmonics matrix [E\ formed as per 
equation (1.1,4). Similarly, for RCWA analyses of TM-polarized incident radiation (131) 
25 according to the formulation of equations (2.3.5) and (2.3,11) or equations (2.3.4) and 

(2.3.10), the required permittivity harmonics e x are the permittivity harmonics s/ 

calculated (410) according to equations (1.1.2) and (1.1.3) and the inverse-permittivity 
harmonics ^-calculated (410) according to equations (2.1.2) and (2.1.3), and the required 
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permittivity harmonics matrices \E\ are the permittivity harmonics matrix [E] formed 
from the permittivity harmonics e x as per equation (1.1.4) and the inverse-permittivity 
harmonics matrix [P] formed from the inverse-permittivity harmonics ^as per equation 
(2.1.4). 

5 As per equations (1.3.7), (2.3.4), (2.3.5) and (2.3.6), the wave-vector matrix 

[A] is dependent on the required, permittivity harmonics matrices \e\ and the matrix [K x ]. 
The matrix [K x ], in addition to being dependent on layer-property parameters (Le. t the 
atmospheric index of refraction no and pitch £), is dependent on incident-radiation 
parameters, Le„ the angle of incidence 9 and the wavelength X of the incident radiation 
10 (131). Therefore, as shown in the flowchart of FIG. 5, according to one embodiment of 
the present invention, ranges 9 m i n to 9 max and ^ m j n to X max , arid increments 66 and 8X are 
determined (617) for the incidence angle 9 and wavelength X, and forwarded from an I/O 
device (805) to the CPU (815). The incident-radiation caching grid {k} is defined as a 
union of two-dimensional coordinates as follows: 

15 {ffHUO?^ +««U mm + 08X) (4.1.3) 

where n and o are integers with value ranges of 

0^n<(6^~a m J/S& , (4.1.4a) 

(4.1.4b) 

(The variable o in equations (4.L3) and (4.1.4b) is not to be confused with the maximum 
20 harmonic order o used in many of the equations above.) Furthermore, the master caching 
grid {fi, K}is defined as a union of coordinates as follows; 

= U K min +iSn r > wot™ +7*o^ D^+kSD, d^+lSd, fi m +mS/3, # min +60, X mm +m8X) 

where ij, k } I, m, n and o satisfy equations (4.1.2a), (4.1.2b), (4.1.2c), (4.1. 2d), (4,1. 4a) 
25 and (4. 1 .4b). Typically, the ranges 9 min to e max and X min to X max are determined (617) 
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based on knowledge and expectations regarding the apparatus (not shown) for generation 
of the incident radiation (131) and the apparatus (not shown) for measurement of the 
diffracted radiation (132). Similarly, the increments 59 and bX are determined (617) 
based on the resolution to which the layer-property parameters n Vy n 0 , A d 9 and p are to be 
determined, and/or the resolution to which the incident-radiation parameters 0 and X can 
be determined. For instance, the increments 8w r , 6no, SD, dd, 5(3, 89, and bX may be 
determined as per the method disclosed in the provisional patent application entitled 
Generation of a Library of Periodic Grating Diffraction Spectra, filed September 15, 2000 
by the same inventors, and incorporated herein by reference. For each point in the master 
caching grid k}, the matrix [A] is calculated (620) by the CPU (815) according to 
equation (1.3.7), (23.4), (2.3.5) or (2.3.6) and cached (425). 

It should be noted that if any of the layer-property parameters n r , n 0i D, d, and 
p, or any of the incident-radiation parameters 9 and X, are known to sufficient accuracy, 
then a single value, rather than a range of values, of the variable may be used, and the 
dimensionality of the master caching grid {\x, k} is effectively reduced. It should also be 
understood that incident-radiation parameter region need not be a hyper-rectangle, and 
the incident-radiation parameter region need not be sampled using a grid. For instance, 
the sampling of the incident-radiation parameter region may be performed using a 
stochastic sampling method. Furthermore, the sampling density of the incident-radiation 
parameter region need not be uniform. For instance, the sampling density may decrease 
near the boundaries of the the incident-radiation parameter region if radiation-incidence 
circumstances near the boundaries are less likely to occur. 

Since the wave-matrix matrix [A] is only dependent on intra-layer parameters 
(index of refraction of the ridges n ft index of refraction of the atmospheric material no, 
pitch D, ridge slab width d, x-offset (3) and incident-radiation parameters (angle of 
incidence 9 of the incident radiation (131), wavelength X of the incident radiation (131)), 
it follows that the eigenvector matrix [W] and the root-eigenvalue matrix [Q] are also 
only dependent on the layer-property parameters n r , n 0 , D 9 d 9 and p, and the incident- 
radiation parameters 9 and X. According to the preferred embodiment of the present 
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invention, the eigenvector matrix [W] and its root-eigenvalue matrix [Q] are calculated 
(647) by the CPU (815) and cached (648) in memory (820) for each point in the master 
caching grid {\x, k}. The calculation (647) of the eigenvector matrices [W] and the root- 
eigenvalue matrices [Q] can be performed by the CPU (815) using a standard 
5 eigensystem solution method, such as singular value decomposition (see Chapter 2 of 
Numerical Recipes, W.H. Press, B.P. Glannery, S .A. Teukolsky and W.T. Vetterling, 
Cambridge University Press, 1986). The matrix [V\, where [F]= [W] [Q], is then 
calculated (457) by the CPU (815) and cached (658) in memory (820). 

The method of use of the pre-computed and cached eigenvector matrices 
10 [^k] ? root-eigenvalue matrices [<2^ K ]> and product matrices [F^ K ] according to the 

present embodiment of the invention is shown in FIG. 6. yse of the cached eigensystem 
matrices [W^ K ], [Q^ 9 and [V^ K ] begins by a determination (505) of the parameters 
describing a discretized ridge profile. In particular, the intra-layer parameters (ie. t index 
of refraction of the ridges n r> the index of refraction of the atmospheric material no, the 
15 pitch D y the ridge slab width d, and the x-offset j$) for each layer, and the incident- 
radiation parameters (Le. } the angle of incidence 9 and the wavelength X of the incident 
radiation) are determined (505) and forwarded via an I/O device (805) to the CPU (815). 
The determination (505) of the discretized ridge profile may be a step in another process, 
such as a process for determining the ridge profile corresponding to a measured 
20 diffraction spectrum produced by a periodic grating. 

Once the intra-layer and incident-radiation parameters are determined (505), 
the cached eigensystem matrices [W^], [<2^ K ], and [V^ K ] for those intra-layer and 
incident-radiation parameters are retrieved (510) from memory (820) for use by the CPU 
(815) in constructing (515) the boundary-matched system matrix equation (1.4.4) or 
25 (2.4.4). The CPU (815) then solves (520) the boundary-matched system matrix equation 
(1 .4.4) or (2.4.4) for the reflectivity R t of each harmonic order from ~o to +o and each 
wavelength X of interest, and forwards the results to an output device (805) such as the 
display (801), printer (803), or the like. 

6. Library Matching Methodology 
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The reflected phase and magnitude signals obtained, in the case of ellipsometry, 
and relative reflectance, in the case of reflectometry, from the profile extraction 
metrology setup are then compared to the library of profile-spectra pairs generated by the 
cached-coupled wave algorithm. The matching algorithms that can be used for this 

5 purpose range from simple least squares approach (linear regression) to a neural network 
based approach that associates features of the signal with the profile through a non-linear 
relationship to a principal component based regression scheme. Explanations of each of 
these methods is explained in numerous excellent text books on these topics such as 
Chapter 14 of "Mathematical Statistics and Data Analysis" by John Rice, Duxbury Press 

10 and Chapter 4 of "Neural Networks for Pattern Recognition" by Christopher Bishop, 
Oxford University Press. 



It should be noted that although the invention has been described in term of a 
method, as per FIGS. 5 and 6, the invention may alternatively be viewed as an apparatus. 

1 5 For instance, the invention may implemented in hardware. In such case, the method 

flowchart of FIG. 5 would be adapted to the description of an apparatus by: replacement 
in step 605 of "Determination of ranges and increments of layer-property variables 
defining array {u} ? and means for determination of maximum harmonic order o" with 
"Means for determination of ranges and increments of layer-property variables defining 

20 array {|u} ? and means for determination of maximum harmonic order o"; the 

replacement in step 617 of "Determination of incident-radiation ranges and increments 
defining array {k}" with "Means for determination of incident-radiation ranges and 
increments defining array {k}"; the replacement in steps 610, 620, 647, and 657 of 
"Calculate ..." with "Means for Calculating ..*"; and the replacement in steps 615, 615% 

25 625,648 and 658 of "Cache ..."with "Cache of 

In the same fashion, the method flowchart of FIG. 6 would be adapted to the 
description of an apparatus by: replacement of "Determination ..." in step 505 with 
"Means for determination . . replacement of "Retrieval ..." in step 5 1 0 with "Means 
for retrieval . . replacement of "Construction ..." in step 515 with "Means for 
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construction . . and replacement of "Solution ..." in step 520 with "Means for solution 

It should also be understood that the present invention is also applicable to off- 
axis or conical incident radiation 1 3 1 (/. e. , the case where <J)^0 and the plane of incidence 

5 140 is not aligned with the direction of periodicity, x , of the grating). The above 
exposition is straightforwardly adapted to the off-axis case since, as can be seen in 
"Rigorous Coupled- Wave Analysis of Planar-Grating Diffraction," M. G. Moharam and 
T. K. Gaylord, I Opt. Soc. Am., vol. 71, 81 1-818, My 1981, the differential equations 
for the electromagnetic fields in each layer have homogeneous solutions with coefficients 

1 0 and factors that are only dependent on intra-layer parameters and incident-radiation 
parameters. As with the case of on-axis incidence, intra-layer calculations are pre- 
calculated and cached. In computing the diffracted reflectivity from a periodic grating, 
cached calculation results for intra-layer parameters corresponding to the layers of the 
periodic grating, and incident-radiation parameters corresponding to the radiation 

15 incident on the periodic grating, are retrieved for use in constructing a boundary-matched 
system matrix equation in a manner analogous to that described above. 

It is also important to understand that, although the present invention has been 
described in terms of its application to the rigorous coupled- wave method of calculating 
the diffraction of radiation, the method of the present invention may be applied to any 

20 diffraction calculation where the system is divided into layers, and where intermediate 
calculation results are only dependent on intra-layer variables. In such case the 
intermediate, intra-layer calculations may be pre-computed and cached. For instance, the 
diffraction calculation may be an approximate method, and/or it may use any of the 
formulations mentioned in the Background of the Invention section, such as integral 

25 formulations, or any other formulations, such as those mentioned in standard texts such 
as Solid State Physics, N.W. Ashcroft and N.D. Mermin, Saunders College, Philadelphia, 
1976, pages 133-134, or Optical Properties of Thin Solid Films, CXS. Heavens, Dover 
Publications, Inc., New York, 1991, or Ellipsometry and Polarized Light, R.M.A. Azzam 
and N.M. Bashara, North-Holland Personal Library, Amsterdam, 1987. Furthermore, the 

30 present invention may be applied to diffraction calculations based on decompositions or 
analyses other than Fourier analysis, such as a decomposition into Bessel functions, 
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Legendre polynomials, wavelets, etc. More generally, the method of the present 
invention may be applied to any diffraction calculation where the system is divided into 
sections, and where intermediate calculation results are only dependent on intra-section 
variables. Again, the intermediate, intra-section calculations may be pre-computed and 
5 cached. For instance, for a two-dimensionally periodic structure, the sections may be a 
regular array of blocks or cubes. Still more generally, the method of the present 
invention may be applied to any calculation where the system is divided into sections, 
and where intermediate calculation results are only dependent on intra-sections variables. 
Again, the intermediate, intra-section calculations may be pre-computed and cached. 

10 The foregoing descriptions of specific embodiments of the present invention 

have been presented for purposes of illustration and description. They are not intended to 
be exhaustive or to limit the invention to the precise forms disclosed, and it should be 
understood that many modifications and variations are possible in light of the above 
teaching. The embodiments were chosen and described in order to best explain the 

15 principles of the invention and its practical application, to thereby enable others skilled in 
the art to best utilize the invention and various embodiments with various modifications 
as are suited to the particular use contemplated. Many other variations are also to be 
considered within the scope of the present invention. For instance: the calculation of the 
present specification is applicable to circumstances involving conductive materials, or 

20 non-conductive materials, or both, and the application of the method of the present 
invention to periodic gratings which include conductive materials is considered to be 
within the scope of the present invention; once the eigenvectors and eigenvalues of a 
wave-vector matrix [A] are calculated and cached, intermediate results, such as the 
permittivity, inverse permittivity, permittivity harmonics, inverse-permittivity harmonics, 

25 pennittivity harmonics matrix, the inverse-permittivity harmonics matrix, and/or the 
wave-vector matrix [A] need not be stored; the compound matrix [V\ 9 which is equal to 
the product of the eigenvector matrix and the root-eigenvalue matrix, may be calculated 
when it is needed, rather than cached; the eigenvectors and eigenvalues of the matrix [A] 
may be calculated using another technique; a range of an intra-layer parameter or an 

30 incident-radiation parameter may consist of only a single value; the grid of regularly- 
spaced layer-property values and/or incident-radiation values for which the matrices, 
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eigenvalues and eigenvectors are cached may be replaced with a grid of irregularly- 
spaced layer-property values and/or incident-radiation values, or a random selection of 
layer-property values and/or incident-radiation values; the boundary-matched system 
equation may be solved for the diffracted reflectivity and/or the diffracted transmittance 

5 using any of a variety of matrix solution techniques; the "ridges" and "troughs" of the 
periodic grating may be ill-defined; a one-dimensionally periodic structure in a layer 
may include more than two materials; the method of the present invention may be 
applied to gratings having two-dimensional periodicity; a two-dimensionally periodic 
structure in a layer may include more than two materials; the method of the present 

1 0 invention may be applied to any polarization which is a superposition of TE and TM 

polarizations; the ridged structure of the periodic grating may be mounted on one or more 
layers of films deposited on the substrate; the method of the'present invention may be 
used for diffractive analysis of lithographic masks or reticles; the method of the present 
invention may be applied to sound incident on a periodic grating; the method of the 

1 5 present invention may be applied to medical imaging techniques using incident sound or 
electromagnetic waves; the method of the present invention may be applied to assist in 
real-time tracking of fabrication processes; the gratings may be made by ruling, blazing 
or etching; the grating may be periodic on a curved surface, such as a spherical surface 
or a cylindrical surface, in which case expansions other than Fourier expansions would be 

20 used; the method of the present invention may be utilized in the field of optical analog 
computing, volume holographic gratings, holographic neural networks, holographic data 
storage, holographic lithography, Zernike's phase contrast method of observation of 
phase changes, the Schlieren method of observation of phase changes, the central dark- 
background method of observation, spatial light modulators, acousto-optic cells, etc. In 

25 summary, it is intended that the scope of the present invention be defined by the claims 
appended hereto and their equivalents. 
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